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The  study  reported  herein  was  conducted  in  1969  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  as  a part  of  the  vehicle 
mobility  research  program  under  DA  Project  1.T061102B52A,  "Research  in 
Military  Aspects  of  Terrestrial  Sciences,"  Task  01,  "Military  Aspects 
of  Off-Road  Mobility,"  under  the  sponsorship  and  guidance  of  the  Re- 
search, Development  and  Engineering  Directorate,  U.  S.  Army  Materiel 
Command . 

The  study  was  conceived  and  carried  out  by- Dr.  Y.  T.  Chou  under  the 
genera]  supervision  of  Messrs.  W.  G.  Shockley  and  S.  J.  Knight,  Chief 
and  Assistant  Chief,  respectively,  of  the  Mobility  and  Environmental 
Division,  and  under  the  direct  supervision  of  Dr.  K.  W.  Wiendieck, 
former  member  of  the  Research  Projects  Group  of  the  Mobility  Research 
Branch.  Dr.  Chou  prepared  this  report. 

, COL  Levi  A.  Brown,  CK,  and  COL  Ernest  D.  Peixotto,  CE,  were  Directors 
of  the  WES  during  the  course  of  this  study  and  preparation  of  this  report. 
Mr.  F.  R.  Brown  was  Technical  Director. 
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Notation 


c Soil  cohesion  determined  from  triaxial  tests 
Cl  Cone  index  (penetration  resistance) 

Second  invariant  of  stress  tensor 
K Yield  stress  in  simple  shear 
P Vertical  load 

R Resistance  of  soil  to  a probe  (constant) 

u,w  Displacements  of  soil  in  r (radial)  and  7.  (vertical)  directions 
respectively 

v Velocity  of  soil  movement 

.> 

V Penetration  velocity 

a,0  Angles  (fig.  1) 

Y Shear-strain  component 

A Interval  of  time 

Strains  in  r (radia.l ) , (circumferential),  and  7,  (vertical) 
directions,  respectively 

0 Normal  soil  stress 

o!.  Deviatoric  stress  tensor 

Oq  Circumferential  stress 

t Shear  stress 

(f>  Angle  of  internal  friction 


Conversion  Factors , Metric  to  British  Units  of  Measurement 


Metric  units  of  measurement  used  in  this  report  can  be  converted  to 
British  units  as  follows: 


To  Obtain 


inches 

square  inches 
inches 

pounds  (force) 
pounds  per  square  inch 
pounds  per  cubic  inch 


square  centimeters 
millimeters 


0.2248 


kilonewtons  per  square  meter 
kilonewtons  per  cubic  meter 


Summary 


The  deformation  and  flow  characteristics  of  a near-saturated  fat 


clay  under  penetration  were  studied.  Penetrations  were  made  in  a 50.8- 
cr.. -diam  mold  with  a circular  cone,  a circular  plate,  and  two  rectangular 
plates  at  speeds  ranging  from  0.004  to  5.6  cm/sec.  The  strain,  strain 
rate,  and  velocity  fields  in  the  soil  were  calculated  from  the  soil  particle 
movements,  determined  by  measuring  the  displacement  of  pellets  embedded 
in  the  soil  before  penetration.  Actual  soil  flow  patterns  determined  from 
velocity  fields  were  studied.  It  was  found  that  if  the  deformation  energy 
of  the  soil  were  assumed  equal  to  the  penetration  energy,  the  former  could 
be  obtained  by  integration  over  the  affected  volume  of  the  deformed  soil. 

The  penetration  resistance  was  thus  computed  on  the.  basis  of  the  Von  Mises 
yield  criterion  and  compared  with  the  measured  penetration  resistance. 
Computed  and  measured  penetration  resistance  values  were  markedly  differ- 
ent; this  casts  some  doubt  on  the  applicability  of  the  Von  Mises  equation 
to  results  of  tests  on  clay  under  penetration  and  on  the  computational 
procedure  employed. 


! 
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EVALUATING  PENETRATION  TESTS  IN  CLAY  FROM  MEASURED 
SOIL  PARTICLE  MOVEMENTS 


Background 


1.  Cone  penetrometer  tests  have  been  used  extensively  at  the  U.  S. 
Army  Engineer  Waterways  Experiment  Station  (WES)  to  evaluate  soil  strength, 
and  empirical  correlations  have  been- established  between  laboratory  cone 
penetration  resistance  and  tire  performance  parameters.  The  main  advantage 
of  the  cone  penetrometer  is  its  inherent  simplicity;  but  it  must  be 
recognized  that  basic  soil  properties,  such  as  cohesion,  internal  friction, 
viscosity,  etc.,  are  not  determined  by  this  instrument.  Therefore,  the 
present  trend  in  penetration  studies  heavily  favors  experimentation  that 
yields  empirical  results. 

2.  Evaluation  of  cone  penetration  resistance  of  soil  by  relations 
based  on  the  theory  of  continuum  mechanics  has  not  been  particularly 
successful.  Most  of  these  analyses  have  been  oriented  toward  static 
bearing  capacity  of  footings  and  have  not  taken  into  consideration  such 
factors  as  viscosity  and  inertia  effects.  Thus,  while  these  theories 
assist  somewhat  in  understanding  the  nature  of  a penetration  test,  they 
do  not  provide  a means  for  evaluating  it. 

3.  Recently,  Yong^  and  Miller^  developed  a new  approach  toward 
rigid  wheel-soil  interaction.  Using  the  concept  of  continuum  mechanics, 
they  approached  the.  problem  by  an  analytical-empirical  method,  visdTo- 

p las ti city , which  had  previously  been  used  successfully  in  analyzing 
metal  processing.  This  method  assumes  that  soil  stresses  under  a wheel 
can  be  calculated  if  the  velocity  field  within  the  soil  is  known.  Yong 
and  Miller  determined  the  velocity  field  by  sophisticated  X-ray  tech- 
niques. Strains  and  strain  rates  were  computed,  and  the  stress  distri- 
bution was  determined  by  a constitutive  equation  (the  Von  Mises  plasticity 
equation)  linking  stresses  and  strain  rates.  Although  this  so-called 
visicpl as ti ci ty  method  has  yielded  no  conclusive  results,  it  seems 


in  contra  :! 


work  of  Yong  and  Miller,  the  three-dimensionality  of  the  problem  is  taken 
fully  into  account.  It  was  hoped  that  the  results  would  shed  some  light 
on  the  fundamental  mechanics  of  penetration  tests  and  allow  the  evaluation 
of  the  potential  of  the  visioplasticity  method  for  more  complex  problems. 


5.  The  purposes  of  this  study  were  to: 

a.  Study  the  deformation  and  flow  characteristics  of  the  soil 
under  penetration, 

b . Examine  the  Von  Mises  equation  in  light  of  the  test  results. 

c.  Evaluate  the  potential  of  the  visioplasticity  method  for 
predicting  soil-vehicle  interaction  theoretically. 

6.  Ten  tests  with  four  different  probes — a circular  cone,  a circular 
plate,  and  two  rectangu! ar  plates — were  conducted  in  a 50. 8-cm-diam*  mold 
of  heavy  clay  at  speeds  varying  from  0.004  to  5.6^ cm/sec.  Only  one  soil 
was  used;  its  moisture  content  averaged  50  percent,  saturation  ranged  from 
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97  to  99  percent,  and  penetration  resistance  ranged  from  38.8  to  129.5  kN/m  . 
The  strain,  strain  rate,  and  velocity  fields  in  the  soil  were  computed  from 
the  movements  of  soil  particles.  The  movements  were  determined  by  measuring 
the  displacement  of  pellets  embedded  in  the  soil  before  penetration. 

7.  Because  of  a lack  of  funds,  the  study  was  shortened  and  simplified 
considerably  witli  respect  to  the  original  plan.  Results  from  a circular 
cone  and  circular  plate  penetrating  the  soil  at  a very  slow  speed 

(0.004  cm/sec)  were,  the  only  data  used  to  compare  measured  and  predicted 
penetration  resistance,  and  results  of  the  other  tests  arc  presented  with- 
out further  analysis. 

Theory tl ca  1_  Background 


Basic  cone  e_p  t 

8.  In  the  basic  concept  of  soil  penetration  problems  in  clay  pic 
seated  in  this  section,  the  validity  of  the  theory  of  perfectly  plastic  solids 
and  the  selection  of  a relatively  hlgh-moisture-content  clay  for  the  test 
program  are  discussed.  Also,  the  morning  of  the  velocity  field  and  the 


*A  table  of  factors  for  converting  metric  to  British  units  of  meac'iverent 
is  given  on  page  ix. 


differences  in  the  stress  conditions  in  cases  of  plane  strain  and  axi- 
symmetry  are  pointed  out,  with  particular  emphasis  on  the  significance 
of  the  circumferential  stress  a. . An  equation  is  formulated  based  on 

U 

the  principle  of  energy  conservation.  This  equation  equates  the  rate  of 
input  penetration  energy  to  the  rate  of  output  soil  deformation  energy, 
and  can  be  used  to  evaluate  the  physical  constants  of  the  material  numer- 
ically or  to  check  available  theories. 

9.  Validity  of  the,  theory  of  perfectly  plastic  solids.  According 

to  the  Mohr-Coulomb  law,  shear  stress  increases  with  normal  stress  on  the 

failure  plane  for  all  but  purely  cohesive  soils.  Drucker  and  Prager 

generalized  the  Mohr-Coulomb  law  to  account  for  all  principal  stresses, 

2 

and  Drucker  demonstrated  that  by  applying  the  normality  rule,  a volume 
increase  must  accompany  any  plastic  strain. 

10.  Upon  application  of  a shear  stress,  soils  undergo  volume  changes 
in  a very  complicated  manner.  Dense  soils  and  overconsolidated  clays  in- 
crease in  volume  upon  shear  stress,  but  normally  consolidated  clays  and 
loose  sands  decrease  in  volume.  Also,  sands  can  undergo  further  plastic 

strain  without  volume  change  once  they  have,  been  strained  far  enough. 

9 

Whitman  pointed  out  that  the  observed  rates  of  expansion  of  dense  sands 
are  far  less  than  those  predicted  from  plasticity  theory. 

11.  In  discussing  the.  importance,  and  necessity  of  the  normality  rule 

9 4 

Whitman  quoted  from  a paper  by  Drucker,  Gibson,  and  Henkel:  'The 

normality  rule  is  much  more  than  simply  an  assumption.  It  is  the  inevi- 
table consequence  of  making  two  less  stringent  and  very  reasonable 
assumptions:  (a)  isotropy  and  (b)  a stable  material.  Also,  the  important 

upper  and  lower  bound  theories  are.  valid  only  so  long  as  the  normality 
ruls  applies.  Without  the  normality  rule,  we  lose  all  theoretical  justi- 
fication for  using  such  results  as  the  Frandt.l-llill  bearing  capacity 
equations . " 

12.  The  general  tendency  in  research  on  plastic,  flow  problems  of 
frictional  materials  is  to  retain  the  normality  rule,  but  seek  for  other 
yield  criteria.  Strain  rates  do  not  remain  normal  to  the  Mohr-Coulomb 
yield  surface,  but  are  normal  to  another  curved  surface.  The  normality 

rule  is  satisfied  if  the  curved  surface  is  the  yield  surface.  Drucker, 

4 

Gibson,  and  Henkel  proposed  that  the  yield  surface  for  soils  should  look 
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like  the  Mohr-Coulomb  surface,  except  that  it  should  be  capped  at  the 
open  end  by  a dome  that  would  expand  and  contract  as  the  volume  of  the 
soil  changed.  In  research  conducted  at  Massachusetts  Institute  of 
Technology  (MIT)  for  the  V7ES,  the  yield  surface  was  assumed  to  be  el- 
liptical and  to  be  able  to  move  along  the  stress  axis  during  plastic 
strain,  the  plastic  volume  changing  in  proportion  to  the  change  in  stress. 

13.  Although  the  concept'  developed  at  MIT  seemed  to  be  reasonable 
and  promising,  conclusive  results  have  not  yet  been  established.  It  is 
felt  that  a correct  plasticity  theory  solving  plastic  flow  problems  for 
frictional  materials  cannot  be  obtained  in  the  near  future.  Since  the 
theory  of  perfectly  plastic  materials  has  been  found  to  be  sound  when 
applied  to  nonfrictionul  materials,  this  study  was  limited  to  soft  clay 
at  a consistency  near  full  saturation,  which  satisfied  the  condition  0=0. 

lU.  The  clay  t-ested  had  a moisture  content  of  5°  percent,  which  was 
only  about  3 percent  short  of  saturation.  The  angle  of  internal  friction 
0 determined  from  triaxial  tests  was  zero.  Although  the  high  moisture 
content  of  the  soil  created  difficulty  and  inconvenience  in  preparing  the 
sample,  it  was  overwhelmingly  advantageous  in  the  theoretical  analysis. 

15 » Velocity  field.  Velocity  fields  in  a soil  mass  show  the  direc- 
tions and  velocities  of  soil  particle  movements',  or  rather  the  flow  or 
rupture  patterns  of  soil.  In  this  study,  the  velocity  fields  were  com- 
puted from  the  measured  soil  particle  movements,  and  the  soil  flow  pat- 
terns under  different  shapes  of  probes  penetrating  at  different  speeds 
were  investigated. 

lu . Stress  conditions.  Stresses  on  an  element  of  soil  in  cylindri- 
cal and  rectangular  coordinates  are  shown  in  figs,  la  and  lb.  In  an 
axisymmetric  case  (fig.  lc),  the  stress  components  are  independent  of  the 
angle  0 , and  all  derivatives  with  respect  to  G vanish.  The  strains 
and  displacements  have  the  relations 
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Cylindrical  coor d in ales 


b . Rectangular  coordinates 


m axisymmetric  case 


1.  Cooruin 


shear-strain  components 


17.  In  the  axi symmetric  case,  displacement  in  the  circumferential 

direction  is  zero,  since  movement  is  confined  to  the  r-z  plane.  But 
strain  in  the  circumferential  direction  is  not  zero  and,  conse- 
quently, stress  in  the  same  direction  exists. 

18.  In  a plane  strain  case,  displacement  in  the  circumferential 
direction  and  strain  e in  the  longitudinal  direction  are  both  zero, 
but  stress  a may  exist.  For  the  plane  strain  case,  (fie.  lb),  equa- 
ti on  1 has  the  form 


19.  In  an  axisyinmetric  case,  it  is  possible  to  determine,  by  the 
condition  of  equilibrium,  the  normal  stress  c and  the  shear  stress 
bn  the  line  inclined  to  the  r 'axis  at  ancle  a in  terms  of  t ,0 
and  the  ancle  a (fig.  1c).  The  equations  can  be  formulated  as 


Stresses  on  an  oblique  line  in  the  r-z  plane,  therefore,  do  not  depend 
on  the  circumferential  stress  o ^ for  the  axisyramctric  case.  It  was 

found  that  equation  3 is  exactly  the  seme  for  a and  t in  the  plane 

. 5 

strain  state  of  stress.  Thus,  the  two-dimensional  Mohr  circle  can  he 


used  to  represent  the  stat 


problems  in  which  the  circumferential  stress 
influence  the  normal  and  shear  stresses  on  the 


maximum  and  minimum  principal  stresses  must  lie  in  the  r-z  plane  for 

the  axially  symmetric  state  of  stress  since  movement  is  confined  to  the 

r-z  plane.  Thus,  0 must  be  the  intermediate  principal  stress.  Accord r 

® 11 

ing  to  the  Harr-Vo'n  Karman  hypothesis,  the  limiting  state  of  stress  at 
yield  is  reached  when  two  of  the  three  principal  stresses—that  is  when 
two  Mohr  circles — are  tangent  to  the  limiting  Mohr  envelope.  Thus,  o0 
must  be  equal  to  either  the  maximum  principal  stress  or  the  minimum 
principal  stress. 

20.  Principle  of  energy  conservation.  Wien  a probe  penetrates  the 
soil  at  constant  speed  V , the  rate  of  input  energy  is  RV  , where  R is 
the  constant  resistance  experienced  by  the  probe.  If  it  is  assumed  that 
the  energy  RV  , which  is  a function  of  response  behavior  of  the  soil,  is 
totally  lost  in  deforming  the  soil,  the  following  energy  balance  equation 
can  be  formulated. 


(Vr  +V0  + Vz  + WrQ  + \z\z  + VW  dV 


The  expressions  shown  at  the  right  are  the  rates  at  which  the  stresses 
deform  the  soil  medium,  and  they  are  integrated  over  the  affected  soil 
volume . • 

21.  In  equation  )(,  the  quantities  R and  V can  be  measured;  and 

components  of  rate  of  strain  tensor  £r  » , and  y^  can  be  computed 

from  the  velocity  field.  With  appropriate  constitutive  equations  that 
link  soil  stresses  to  strain  rates,  the  components  of  stress  tensor  o > 
°0  ’ T0-  * etc*>  can  t>e  obtained  in  terms  of  strain  rate,  and  the  inte- 
gration can  be  carried  out.  For  tests  at  very  slow  penetration  speeds, 
the  Von  Mises  equations  for  perfectly  plastic  solids,  were  used  as  the 
constitutive  equation.  The  results  of  these  tests  are  presented  for 
future  use,  and  the  velocity  fields  are  plotted. 

Equations 

22.  Computation  of  strain  rate  fields.  For  an  axi symmetric  case, 
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rates  of  strain  in  r , 6 , and  z directions 
respectively 

velocities  in  r and  z directions,  respecti 


rate  of  shear-strain  in  r-z  direction 


For  a plane  strain  case 


23.  Equations  5 and  C are  valid  only  for  small  deformations  and  are 
in  a strict  sense  not  accurate  enough  for  problems  with  large  strains. 
However,  in  view  of  present  crude  measuring  techniques  and  other  idealiza- 
tion of  soil  properties,  it  appears  unrealisitc  to  include  nonlinear  terms 
in  equations  5 and  6 to  account  for  larger  strains.  Nevertheless,  these 
equations  were  used  for  computing  strain  rates. 

Lcvy-Von  Hiscs  plasticity  equations.  The  Levy-Von  Miser  equa- 
tions are: 


stress  deviation 
respec tively 


It  is  frequently  convenient  to  decompose  the  stress  tensor  into  a s 
cal  part  corresponding  to  the  roan  normal  stress  and  a clevintoric 
The  stress  deviation  is  defined  as  the  tensor  with  the  normal  cor.? -on 
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25.  These  equations  assume  that  (a)  the  rate  of  strain  is  propor- 
tional to  the  stress  deviation,  (b)  the  material  is  incompressible  and 
nonviscous,  and  (c)  the  condition  - K is  never  violated ‘during  flow. 
(J2  is  the  second  invariant  of  the  stress^tensor  and  K is  the  yield  ' 
stress  in  simple  shear. ) This  last  condition  implies  that  work-hardening 
phenomena  do  not  exist. 

2 6.  Since  the  test  soil  was  near  saturation  and  could  be  assumed  to 
be  incompressible,  hydrostatic  stresses  would  not  cause  any  soil  straining 
or  energy  dissipation,  i.e.  the  hydrostatic  stress  did  not  influence  the 
deviatoric  stress  behavior.  Hence,  in  energy  computations,  as  in  equa- 
tion U,  it  was  necessary  to  consider  only  the  deviatoric  component  of 
stress,  which  can  be  represented  by  a!.  , the  deviatoric  stress  tensor. 
Material  constants 


2'f.  In  the  theoretical  analysis  of  soil  mechanics  problems,  it  is 
very  difficult  to  evaluate  the  material  constants  from  laboratory  test 
data.  These  constants  depend  on  a number  of  factors , probably  the  most 
important  being  the  difference  in  the  stress  systems,  i.e.  the  laboratory- 
determined  material  constants  are  based  on  ideal  loading  conditions  that 
usually  do  not  represent  the  real  conditions  under  a complex  state  of 
stress.  If  it  is  assumed  that  the  laboratory  triaxial  test  results  can  be 
used  to  determine  soil  strength,  K = (ci^  - a )/»^3  when  the  Von  Mises 
yield  condition  is  used.  K ~ v^-en  Tresca  yield  condition 

is  used. 

28.  In  this  study,  the  Von  Mises  yield  criterion  was  used.  The  K 
value  was  not  determined  directly  from  the  triaxial  tests,  which  are  dif- 
ficult and  expensive  to  run  at  such  high  moisture  contents,  but  by  extrap- 
olations from  previous  triaxial  and  cone  penetration  test  results  with  the 


which 


c.  Compute  the  K value  from  K = 2c  //3  . 

2 

The  penetration  resistance  of  the  test  soil  was  found  to  be  55  kll/m  , 
yielding  a cohesion  of  U . U klf/iri  . K was  computed  as  5.1  kll/in  . 

29 • Based  on  results  of  numerous  tests  on  both  fat  and  lean  clays 
with  moisture  contents  of  U5  percent  or  less  (no  tests  were  conducted 
at  50  percent  moisture  content  comparable  to  the  present  study),  the 
relation  shown  in  equation  8 was  formulated.  Although  the  error  caused 
by  the  linear  extrapolations  beyond  the  test  range  is  not  exactly  known 
the  K value  thus  obtained  should  be  fairly  reasonable.  It  is  very 
doubtful,  however,  if  this  K value  is  exactly  the  yield  stress  that 
occurred  in  the  soil  under  penetration. 


Laboratory  Test 


Soil  and  it 


30.  The  fat  clay  used  in  this  study  was  river  alluvium  obtained 
near  Vicksburg,  Miss.,  and  was  classified  CH  according  to  the  Unified 
Soil  Classification  System.  Gradation  and  classification  data  are  shown 
in  fig.  2. 

31.  The  soil  was  air  dried  to  a uniform  moisture  content,  crushed 
with  a hammer  mill,  and  sieved  through  a 0.3-cin  screen  before  being  mixed 
with  water.  The  soil  and  water  were  thoroughly  blended  in  a pug  mill 
until  a uniforn  consistency  was  obtained,  after  which  it  was  stored  in  an 
airtight  container  for  about  2h  hr  before  being  molded  into  a sample. 
Sample  preparat i on 

32.  A test  sample  was  prepared  in  a 50. 8- cm- di am  mold  that  could  be 
separated  vertically  into  two  halves.  Processed  soil  was  compacted  in  the 
mold  with  a mechanically  operated,  327-N,  wedge-shaped  drop  hammer  falling 
15.2  cm.  A l-cm- thick  rubber  mat  was  placed  over  the  surface  of  the  soil. 


but  it  was  lifted  occasionally  to  prevent  its  sticking  to  the  soil  and 
blocking  the  free  release  of  air.  A sample  was  made  up  of  20  layers,  each 
receiving  1)00  blows  of  the  compaction  hammer  uniformly  over  the  surface. 
With  this  procedure,  it  is  believed  that  the  sample  achieved,  uniformity, 
and  a maximum  amount  of  air  was  expelled  from  the  soil , 


Hi  u.w«ii  ii'pjn'i  Bi[pwi^iiiriTTTwr^niiiTrrTrniriir|i|>lwrTIWITI>1l1ITIinnTM*FM*M!*l"M‘**w 


Testing  procedure 

33.  After  the  soil  had  been  compacted,  the  two  halves  of  the  mold 
were  slightly  disengaged,  and  the  sample  was  cut  vertically  in  half  with 
a piano  wire.  The  mold  was  disassembled  and  laid  on  a table  with  the  flat 
soil  surface  up. 

3^.  Straight,  orthogonal  grid  lines  were  drawn  at  1, 27—cm  intervals 
on  the  soil  surface  with  a sharp  pencil.  At  locations  where  the  grid 
lines  would  be  destroyed  by  penetration,  small  colored  pins  were  inserted. 
The  complete  pattern  of  the  grid  lines  then  was  traced  on  a transparent 
plastic  sheet.  After  the  soil  surfaces  were  oiled  lightly,  the  two  halves 
of  the  mold  were  put  back  together,  and  the  penetration  tests  were  made 
near  the  center  of  the  mold. 

35.  After  penetration,  the  mold  was  disassembled  and  the  halves  were 
separated  again.  The  deformation  pattern  of  the  grid  lines  was  traced  on 
the  same  plastic  sheet.*  Where  the  lines  were  destroyed  by  penetration, 
the  areas  could  be  reconstructed  readily  by  the  new  positions  of  the 
colored  pins. 

Accuracy  of  test  method 

36.  The  test  method  used  in  this  study  was  rather  crude.  Factors 
that  affected  the  reproducibility  of  test  results  were  found  to  be: 

a.  Variation  (within  +2  percent)  in  moisture  content  of  the 


samples. 

b.  Nonuniformity  of  the  samples. 

c.  Variation  in  the  penetration  technique  (i.e,  a slight  eccen- 


tricity of  the  penetration  could  not  always  be  avoided). 


In  two  different  soil  samples  prepared  and  penetrated  under  identical  con- 
ditions, the  penetration  resistance  values  were  very  close,  but  the  deformed 
grid  lines  (or  the  soil  flow  patterns)  were  not  exactly  the  same.  In  other 
words,  the  displacement  of  a soil  particle  at  a particular  location  in  one 
soil  sample  was  not  the  same  as  in  the  same  location  in  a second  sample. 
However,  this  did  not  seen  to  constitute  a problem  in  this  study,  because 
the  computations  were  based  on  integration  over  the  total  volume  of  the 


* In  both  axisynnetric  and  plane  strain  cases,  soil  particles  moved  only 

in  the  cut  plane. 
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deformed  soil  and  did  not  depend  on  individual  particle  movements. 


Presentation  of  Test  Results  and  Their  Analysis 


37*  The  results  of  this  study  are  presented  in  three  parts:  (a)  lab- 

oratory test  results  that  show  the  flow  characteristics  and  particle  move- 
ment of  the  soil  under  penetrations,  (b)  the  computer  program,  and  (c)  the 
theoretical  prediction  of  penetration  resistance.  The  penetration  resist- 
ance was  predicted  by  assuming  the  deformation  energy  of  the  soil  to  be 
equal  to  the  penetration  energy,  and  an  integration  of  the  deformation 
energy  over  the  total  affected  volume  of  the  deformed  soil  was  carried  out 
to  obtain  the  penetration  resistance.  The  computed  penetration  resistance 
then  was  compared  to  the  measured  cone  penetration  resistance. 

Laboratory  test  results- 


38.  Representative  laboratory  test  results  are  shown  in  table  1 and 
plates  1-10.  IJonhomogeneity  of  the  test  soils  could  not  always  be  avoided, 
so  the  deformation  patterns  on  opposite  sides  of  the  probe  were  usually 
dissimilar.  This  was  reasonable  because  the  load  tended  to  deform  most  on 
the  side  that  offered  the  least  resistance.  Occasionally,  the  magnitude 
of  displacement  at  a distant  point  was  found  to  be  greater  than  that  at 
a point  closer  to  the  probe. 

' 39 • Test  results  for  a 30-de’g  circular  cone  with  a base  area  of 

2.3  cm  cm  in  diameter)  are  presented  in  plates  1—4 , Penetration 

speed  was  5*6  cm/sec,  the  highest  used.  (The  pins  near  the  cone,  in 
plate  1 were  placed  after  the  test  to  provide  a better  contrast  in  the 
photograph;  this  was  found  later  to  be  unnecessary.)  The  affected  zone 
was  rather  small,  and  soil  beneath  the  tip  of  the  cone  seems  not  even  to 
have  moved.  The  stream  lines  (flow  path")  of  soil  particles  are  shown 
in  plate  2.  Plate  3 is  the  complete  trace- of  soil  particle  movements  at 
the  grid  points  from  the  time  when  the  cone  was  Just  touching  the  soil 
surface  until  the  end  of  the  test.*  The  displacement  vectors  of  soil 
of  soil  particles  at  grid  points  around  the  cone  are  presented  in  plate  1;; 


* This  result  v r.  obtained  by  as  a um.pt  5.  on  of  steady  state,  which  is  di 
cussed  in  paragraph  . 
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the  starting  and  ending  points  of  each  arrow  show  the  original  and  final 
positions  of  the  soil  particles. 

HO.  Results  of  tests  with  a 5.1-cm-diam  circular  plate  moving  at  an 
extremely  slow  speed  (O.OOH  cm/sec)  are  presented  in  plates  5-7 , and  with 
a 3.8-  by  22.8-cm  rectangular  plate  at  the  same  speed  in  plates  8-10. 

(The  pins  shown  in  plates  5 and  8 vere  placed  before  the  tests.)  Under 
this  slow  penetration  speed,  the  soil  was  generally  assumed  to  behave  as  a 
perfectly  plastic  material,  an  assumption  that  did  not  necessarily  apply 
to  fast  penetration  speeds. 

Hi.  The  portion  of  soil  that  moved  with  each  plate  is  termed  the 
soil  nose  in  this  report.  Soil  noses  under  the. circular  plate  penetrat- 
ing at  four  different  speeds  are  shown  in  figs.  3b-3e,  and  under  the 
5.1-  by  25.H-cm  rectangular  plate  at  two  different  speeds  in  figs,  hb 
and  He . Figures  3a  .and  Ha  show  the  locations  of  pellets  embedded  in  the 
soil  and  the  original  grid  lines  before  penetration.  The  soil  nose  formed 
under  the  5.1-  by  25.H-cn  rectangular  plate  by  penetration  at  a very  slow 
penetration  speed  (fic.  Hb)  had  a triangular  shape,  but  changed  to  cir- 
cular as  the  speed  was  increased  to  the  standard  3.1  cn/sec  (fig.  He). 

H2.  Soils  in  the  vicinity  of  the  soil  noses  experienced  large  shear 
strains.  Because  of  the  adhesion  force  between  the  plates  and  the  soil, 
the  portion  of  soil  near  the  plates  did  not  move  away,  but  moved  vertically 
downward  with  the  plate  like  a rigid  body.  This  may  be  verified  by  the 
position  (relative  to  the  plate)  of  pellet  IIo.  H,  which  was  1.27  cm  below 
the  center  of  the  plate  before  and  after  penetration  in  all  tests.  Pellet 
No.  1 did  not  move  away  in  any  of  the  tests;  pellet  IIo.  2 moved  slightly 
in  some  tests  in  which  the  soi.1  conditions  on  the  right,  side  of  the  mold 
vere  apparently  weaker  than  on  the  left.  Pellet  No.  8,  which  was  placed 
at  a depth  half  the  width  of  the  plate  (2.5H  cm)  below  the  center  of  the 
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of  the  mold  were  heaved  slightly,  indicating  the  mold  actually  vas  not 
large  enough  for  the  deformation  patterns  to  be  complete.  The  soil 
evidently  was  failed  by  shear  stresses,  the  failure  lines  showing  clearly 
in'  plate  8.  Although  the  soil  was  failed  by  shear  stresses  in  the  axisym- 
metric  case  (plates  5 and  7 ) , no  rupture  lines  can  be  seen?  -This  soil 
responded  to  penetration  as  viscous  fluids  might,  and  the  soil  deforma- 
tions were  limited  to  small  areas  near  the  plate.  Despite  these  differ- 
ences, the  unit  penetration  resistances  for  these  cases  were  not  very  dif- 
ferent from  each  other,  as  shown  in  table  1.  When  the  flow' patterns  in 
plates  1 and  5 are  compared,  it  can  be  seen  that  penetration  speed  had  no 
effect  on  soil  deformation  patterns  in  the  axisymmetric  case.  In  the  plane 
strain  case,  however,  the  failure  lines  are  still  seen  at  a penetration 
speed  of  3.1  cm/sec,  but  not  as  clearly  as  at  slower  penetration  speeds. 

Mt.  The  soil  f-low  patterns  in  the  zone  near  the  penetration  hole 
were  not  exact  because  the  soil  here  tended  to  flow  toward  the  hole  by 
gravity.  This  lateral  soil  movement  occurring  during  and  after  the  pene- 
tration gave  false  information  for  computed  results  of  strain  rate  fields. 
The  magnitude  of  this  error,  however,  could  hardly  be  evaluated. 

1*5.  To  trace  the  continuous  movements  of  soil  particles  under  pene- 
tration, the  deformed  grid  lines  are  recorded  as  a function  of  time.  This 
can  be  done  by:  (a)  the  X-ray  method  used  by  Yong^  and  Miller,^  and 

(b)  the  measurement  of  deformed  grid  patterns  in  identical  soil  samples 
penetrated  at  different  depths.  The  first  method  could  not  be  used  for 
this  study  because  X-ray  facilities  for  testing  large  samples  were  not 
available  at  the  WES . The  second  method  vas  not  used  because  a large 
number  of  tests  are  required,  and  identical  soil  samples  are  very  diffi- 
cult to  prepare.  Instead,  steady  state  of  movement  of  soil  particles  vas 
assumed;  this  allowed  the  continuous  soil  particle  movement  to  be  deduced 
from  the  final  positions  of  the  deformed  Grid  lines.  Actually  the  valid- 
ity of  the  use  of  the  superposition  technique  to  trace  the  entire  history 
of  soil  particle  movements  is  tires  assumed. 

1*6.  The  grid  lines  in  a test  sample  being  penetrated  by  a.  cone  are 
shown  schenaticaly  in  fig.  5*  The  horizontal  lines  before  penetration 
are  shown  in  fig.  5a,  and  the  deformed  patterns  at  successive  times  in 
figs.  51>,  5c,  and  5d.  Under  assumed  steady  state,  line  3-d  should  have 
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the  same  shape  as  lines  2-c  and  1-b;  also,  line  2-d  should  have  the  same 
shape  as  line  1-c.  As  the  cone  penetrates  further  into  the  soil,  the 
deformed  lines  surrounding  the  cone  should  have  the  same  shape  as  those  • 
at  At  time  ago.  In  other  words,  the  patterns  of  deformed  lines  are 
index>endent  of  depth,  just  as  the  waves  surrounding  a sailboat  do  not 
change  as  long  as  the  boat  is  moving  at  constant  speed. 

47.  Experimental  results  obtained  at  the  WES°  have  shown  that 
penetration  resistance  of  clay  does  not  depend  upon  penetration  depth, 
except  for  small  variation  due  to  overburden.  The  strength  of  the  soil 
surrounding  the  probe  does  not  change  as  the  probe  penetrates  the  soil 
so  the  soil  can  be  assumed  to  be  a uniform  and  homogeneous  medium  with 
constant  strength.  Hence,  the  assumption  of  steady  state  for  the  cases 
investigated  herein  seems  to  be  valid. 

Computer  program  and  results 

48.  A conqiuter  program  (Appendix  A)  was  prepared  to  be  used  for  pre- 
dicting penetration  resistance  of  the  soil.  It  vas  written  in  three  parts: 

a.  Computation  of  velocity  field. 

b.  Computation  of  strain  rates  at  grid  points. 

c.  Computation  by  the  Von  Mises  equation  of  the  stress  com- 
ponents at  each  grid  point  and  the  rate  of  stresses  deforming 
the  soil 


Based  on  the  energy  balance  equation  4,  the  validity  of  the  Von  liises  equa- 
tion was  checked.  The  numberical  value  of  soi 1 yield  stress  was  deter- 
mined by  laboratory  tri axial  tests  previously  conducted. 

49.  The  essential  work  in  the  program  was  to  compute  th  .■  strain  rates 
at  the  grid  points.  Equations  5 and  7 were  used  for  this  purpose  with  the 
input  data  measured  from  test  results. 

50.  The  original  grid  lines  (solid)  and  the  stveam  lines  (dashed) 
are  shown  in  fig.  6 a;  fig.  6b  shows  the  cample  trace  of  soil  particle 
movement  at  each  grid  point,  from  time  t - C to  t ~ JAt  where  At 

is  the  time  required  for  the  cone  to  move  through  the  distance  D . The 
displacements  in  r and  7 directions,  as  a function  of  time  t , are 
presented  in  figs,  6c  and  6d , respectively.  There  are  various  ways  to 
compute  the  velocity  field  and  strain  rate  field.  That  used  to  compute 
a velocity  field  is  as  follows: 


a.  The  locations  of  deformed  grid  points,  i.e.  a'  , b'  , c'  , 
d'  , e'  , and  f'  (fig.  6a),  were  first  read  into  the  com- 
puter. The  measurements  were  taken  with  reference  to  a 
fixed  coordinate  system. 

b.  The  displacements  u and  w were  computed  according  vo  the 

definition  (u  or  w)^  = total  displacement  up  to  tine  t 

minus  displacement  at  tine  t . 

c_.  An  n-th  order  polynomial  based  on  the  least  square  tech- 
nique was  generated  by'  the  computer  and  represented  the 
displacement  versus  time  curves.  The  slopes  of  the  curves 
at  time  t ~ At  , 2At  , 3At  , ^At. , and  5At  resulted  in 
velocities  cc: responding  to  points  e*  , d'  , c*  , b'  and 
a*  , respectively  (fig.  6a).  It  was  found  that  polynomials 
with  an  order  of  8 would  fit  most  data  with  an  acceptable 
degree  of  accuracy. 


To  compute  a strain  rate  field,  the  following  steps  were  taken: 

a.  Based  on  values  of  velocities  at  points  a'  ,b'  ,c'  ,d'  , 
and  c'  along  a stream  line,  velocities  at  points  b"  , 

c"  , d"  , and  e"  were  obtained.  These  points  were  actu- 
ally the  intercepts  of  the  stream  line  on  the  unde formed 
horizontal  grid  lines. 

b.  After  these  velocities  had  been  obtained  for  all  stream 

lines , the  velocities  at  each  intersect! on  on  the  undeformod 
grid  were  generated  by  the  computer,  and  the  horizontal  and 
vertical  components  of  the  velocities,  as  well  as  the  rate 
of  change  of  velocities  in  both  directions,  were  determined 
by  the  least  square  method  in  conjunction  vith  polynomials . 
Thus,  the  strain  rate  tensor  l.,  at  each  grid  point  was 
obtained.  1 1 


52.  Difficulty  was  encountered,  however,  in  using  the  least  square 
procedure  in  paragraph  5^-h.*  According  to  computer  print-outs,  velocities 
did  not  always  increase  or  decrease  nonot on i c ally , but  fluctuated  in  one 
direction.  For  example,  velocities  along  a vertical  grid  line  below  the 
probe  did  not  always  decrease  with  depth.  It  was  found  that  the  varia- 
tion in  the  order  of  polynomials  was  very  sensitive  to  the  variation  in 
velocity  values.  The  correct  order  for  a given  set  of  data  could  not  be 
determined  analytically  by  any  method  except  trying  a series  cf  different 
values.  In  this  program , orders  from  3 to  lb  (less  than  the  data  points) 
were  tried,  and  the  one  that  gave:  the  best  approximation  was  selected  in 
the  computation. 

53.  A velocity  field  shove  the  flow  pattern  of  soi]  under 


penetration.  The  velocity  fields  beneath  probes  at  the  moment  penetra- 
tions were  completed  are  depicted  in  plates  11-16.  (Since  steady-state 
conditions  were  assumed,  the  same  velocity  field  should  be  observed  at  any 
other  stage  of  penetration.)  The  arrows  shown  indicate  the  magnitude  of 
velocity  and  direction  of  movement  of  soil  particles,  and  the  dashed  lines 
show  the  stream  lines.  These  lines  indicate  the  directions  of  maximum 
shear  stress.  Along  the  stream  lines,  velocities  of  soil  particles  were 
relatively  large  close  to  the  probe,  but  became  smaller  as  the  distance 
increased  from  the  probe.  In  other  words,  the  velocity  decreased  along 
each  stream  line.  The  rates  of  change  of  velocity  (or  strain  rate)  along 
the  stream  lines  at  slow  penetration  speed  (plates  11,  13,  and  15)  were 
observed  to  be  smaller  than  those  at  faster  penetration  speeds  (plates  12, 
lH,  and  16).  Since  the  deviatoric  stress  is  proportional  to  the  strain 
rate,  according  to  the  Von  Mises  concept,  the  stresses  along  equivalent 
stream  lines  should  be  greater  under  faster  penetrations  than  under  slower 
ones.  This  is  reasonable  because  experimental  results  revealed  that  the 
penetration  resistance  increased  with  increasing  penetration  speed. 

5^.  It  is  interesting  to  note  that  soils  under  penetration  do  not 
move  along  a unique  slip  surface  in  the  velocity  fields  as  most  theoretical 
analyses  have  assumed  (fig.  7 ) • It  usually  is  assumed  that  at  the  time  of 


Soil  movements  beneath  a footing 


the  case,  however,  for  the  cases  in  this  study:  The  velocity  varied  along 

the  stream  line  and  also  among  stream  lines;  also,  the  grid  squares  were 
distorted  after  movement  (plates  1,  5»  and  8)f  indicating  that  there  were 
relative  displacements  among  soil  particles  that  moved  along  infinite 
numbers  of  surfaces  rather  than  along  a unique  surface. 

55.  It  is  important  to  realize  that  the  assumption  of  rigid  body 
movement  of  a soil  mass  beneath  a footing  applies  to  the  condition  of 
incipient  failure  or  limit  equilibrium,  i.e.  the  initial  movement  of  the 
soil  mass  at  the  beginning  moment  of  failure.  With  continued  movement, 
as  in  this  study,  the  soil  mass  may  deviate  from  the  simple  rigid  body. 

56.  The  computed  velocity  field  beneath  a 3.8-  by  22.8-cm  rectangu- 
lar plate  moving  downward  at  a very  slow  speed  (O.OOl)  cm/sec)  (plate  15) 
can  be  compared  with  the  rupture  lines  in  the  soil  sample  after  the  test 
(plate  8).  It  appears  that  the  computed  directions  of  soil  particle  move- 
ment checked  quite  well  with  the  observed  direction  or  rupture  lines,  indi- 

> 

eating  that  the  computation  method  for  velocity  fields  employed  in  this 
study  was  reasonably  accurate. 

57.  In  the  velocity  field  beneath  a 5.1-  by  25Ji-cm  rectangular 
plate  moving  downward  at  0.  speed  of  3.1  cm/sec  (plate  l6),  the  stream 
lines  (velocity  field)  appear  to  be  different  from  those  in  plate  15.  At 
slow  penetration  speed  (plate  15) » the  soils  moved  away  from  the  probe  in 
a markedly  horizontal  direction;  whereas  at  faster  speed  (plate  l6) , soil 
near  the  probe  moved  in  a markedly  vertical  direction. 

58.  Velocity  fields  beneath  a 5 . 1-cm-dian  circular  plate  at  pene- 
trating sj>eeds  of  0.004  and  5*6  cm/sec  are  shown  in  plates  13  and  lH, 
respectively.  The  velocity  fields  differed  under  these  two  speeds > os  did 
those  in  plates  15  and  16,  and  the  same  discussions  apply . As  the  pene- 
tration speed  increased,  the  stream  lines  generally  became  steeper,  and 
the  distance  from  the  probe  to  the  ends  of  the  flow  lines  decreased. 

59.  Plates  11  and  12  are  similar  to  plates  13  end  3)t , except  that  a 
5. )+-cr.i  circular  cone  war.  used.  Because  of  the  sharp  apex  angle  of  the 
cone,  soils  beneath  the  cone  tip  were  hardly  moved.  Otherwise,  the  dis- 
cussions previously  mode  for  other  types  of  probe  are  also  applicable  to 
results  with  this  one. 


Theoretical  prediction 
of  penetration  resistance 

60.  The  procedure  for  computing  the  penetration  resistance  based  on 
the  measured  strain  rate  fields  is  given  in  this  section.  The  computer 
results  were  limited  to  the  axi symmetric  coses  at  very  slow  penetration 
speeds  (Nos.  1 and  2 in  table  l). 

61.  With  the  velocity  components  at  each  grid  point  known,  the  com- 
ponents of  strain  rate  can  be  computed;  with  the  proper  constitutive 
equation  for  the  soil,  the  total  resistance  of  the  soil  can  be  evaluated 
by  equation  k (the  energy  balance  equation)  and  theoretically  should  be 
equal  to  the  measured  penetration  resistance.  In  using  equation  1+  to 
carry  out  the  integration,  only  the  affected  soil  volume  beneath  the  probe 
was  considered;  the  effect  of  soil  above  the  probe  was  discarded.  This 
simplification  was  assumed  as  a direct  consequence  of  the  steady-state 
assumption ; it  reduced  the  effect  of  the  soil  above  the  probe  to  simple 
overburden  pressure,  which  was  insignificant  in  the  investigated  cases 
since  the  penetration-  resistance  in  clay  remained  almost  constant  after 
full  penetration  of  the  probe.  Also,  the  soil  above  the  probe  tended  to 
move  toward  the  penetration  hole,  resulting  in  unrealistic  and  false 
velocity  and  strain  rate  fields  insofar  as  the  prediction  of  penetration 
resistance  was  concerned. 

62.  To  save  computer  time,  it  was  assumed  that  the  deformation 
energy  in  the  soil  under  penetration  was  symmetrical  with  respect  to  the 
probe.  Consequently,  the  computations  were  carried  out  for  one  half  of 
the  mold  only.  In  the  analysis,  the  Levy-Von  Mines  plasticity  equation 
was  used.  The  soil  yield  stress  K was  taken  as  5«1  kli/n'  (paragraph  20) 
The  soil  resistances  computed  from  the  measured  strain  rate  fields  were 
very  much  larger  than  the  measured  penetration  resistances,  as  shown  in 
the  following  tabulation. 


ance 


63.  If  the  theory  of  plasticity  is  assumed  to  be  valid  for  friction- 
less materials,  and  if  the  measuring  technique  and  the  computed  strain 
rates  are  assumed  to  be  reasonably  accurate,  questions  may  arise  as  to 
whether  the  numberical  value  of  soil  yield  stress  K;  used  in  computations 
is  properly  chosen.  In  the  Levy-Von  Mises  equation,  K is  defined  as  the 
yield  stress  in  shear,  or  1//3  yield  stress  in  tension  or  compression.  It 
is,  of  course,  questionable  whether  the  yield  stress  determined  from  the 
laboratory  triaxial  test  can  represent  the  yield  stress  under  the  complex 
state  of  stress  in  the  soil  sample  under  penetration.  If  the  computation 
including  the  simplifying  assumptions  is  reasonable,  the  real  soil  yield 
stress  K under  penetration  may  be  computed  by  equation  h.  Substituting 
the  Levy-Von  Mises  plasticity  equation  6 into  equation  *4  yields 


Since  R and  V are  measured  and  the  components  of  the  strain  rate  tern 
sor  can  be  computed,  the  soil  yield  stress  K under  penetration  can  be 
computed  from  equation  9-  These  values  are  listed  in  the  last  column  of 
the  tabulation  in  paragraph  62. 

61).  The  large  discrepancies  between  the  measured  and  computed  pene- 
tration resistance  values  in  these  two  tests  may  be  explained  by : / 

a.  The  soil  -was  nonuniform,  i.e.  the  deformation  patterns  at 
the  cut  surface  did  not  veil  represent  those  on  other 
surfaces. 


The  Levy-Von  Mises  equation  was  inadequate  for  describing 
the  soil  behavior  under  penetration,  i.e  according  to  the 
Von  Mises  equation,  it  is  reasonable  io  believe  that  the 
rate  of  strain  is  proportional  to  the  stress  deviation  dur- 
ing flow , but  it  is  questionable  whether  the  yield  condition 
2 

J„  e K is  always  valid  during  flow.  Also,  the  soil  sample 
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derivatives  of  the  polynomial  would  deviate  very  much  from 
those  of  the  displacement-time  relation  due  to  the  oscillat- 
ing nature  of  the  high-order  polynomial.  This  may  best  be 
explained  by  an  example.  The  solid  line  connected  by  crosses 
in  fig.  8a  shows  the  measured  displacement-time  relation  of 
a soil  particle  in  the  vertical  direction.  A seventh-order 
polynomial  was  generated  based  on  these  data  points  and  was 
plotted  in  the  line  connected  with  circles.  The  portion  of 
the  curve  between  tine  0 to  time  5At  is  plotted  at  an  en- 
larged scale  in  fig.  8b.  The  horizontal  line  was  actually 
approximated  by  an  oscillating  curve  with  the  difference 
less  than  0.023  mm.  When  the  derivatives  of  the  polynomial 
were  tshen  at  various  times,  however,  the  values  were  not 
zero.  As  a result,  nonzero  velocities  or  nonzero  strain 
rates  existed  in  this  time  interval.  In  other  words,  posi- 
tive deformation  energy  was  computed  in  the  portion  of  the 
soil  in  which  measured  soil  movements  were  actually  zero. 

The  deformation  energy  was  always  positive  because  the  com- 
ponents of  the  strain  rate  tensor  in  equation  9 were  squared 
before  summing  for  every  small  soil  volume-  In  some  cases, 
when  the  noil  particles  were  undergoing  rigid-body  movements, 
i.e.  the  soil  particles  moved  at  constant  velocity  without 
interval  energy  dissipation,  positive  deformation  energy 
would  be  computed  for  the  same  reason  described  above. 


Among  all  these  reasons,  it  is  believed  that  the  errors  caused  by  high- 
order  polynomial  approximations  were  most  serious.  To  overcome  this  dif- 
ficulty, the  computer  program  can  be  modified  as  described  below. 

65.  The  measured  displacements  should  be  smoothed  before  being  fed 
into  the  computer.  ( The  dashed 'line  in  fig.  8a  is  a smoothed  displacement' 
time  curve.)  In  the  time  interval  between  zero  and  5At  , the  velocities 
would  be  assumed  equal  to  zero.  A polynomial  with  very  low  order  or  a 
simple  function  would  be  accurate  enough  to  approximate  this  displacement- 
time  curve  without  oscillation.  It  should  be  pointed  out  that  the  soil 
displacemnts  were  actually  very  small,  the  differences  between  the  smooth 


curve 


less  than  1 min  were  determined  quite  arbitre 
believe  that  the*  smoothed  curve  would  better 


•e  con srda red  a continuous  medium 


Conclusions 


Based  on  the  test  data  and  computer  results,  it  is  concluded  that: 

a The  soil  patterns  under  penetration  can  be  drawn  based  on  the 
measured  soil  particle  movements  (paragraphs  50-59). 

b.  The  soil  flow  patterns  in  two-dimensional,  plane  strain  cases 
are  different  from  those  in  axi symmetric  cases.  In  both 
cases,  the  shape  of  the  stream  (or  flow)  lines  change  with 
change  in  penetration  speeds.  In  general,  at  the  high  pene- 
tration speeds  of  this  study,  the  rate  of  change  in  the  slope 
of  the  stream  lines  is  greater  and  the  distances  from  the 
ends  of  the  stream  lines  to  the  probe  become  smaller  than  at 
slow  speeds.  Also,  at  the  high  penetration  speeds,  the  rate 
of  velocity  change  (or  strain  rate)  along  stream  lines  is 
greaterj  stresses  in  the  soil  are  consequently  greater  to 
counterbalance  the  higher  penetration  force  (paragraph  53). 

c_.  The  clay  beneath  the  moving  probe  does  not  move  entirely  as 
a rigid  body  along  a unique  failure  surface;  rather  it  moves 
along  infinite  numbers  of  failure  surfaces  (paragraph  5*0. 

d_.  The  penetration  resistances  computed  from  measured  soil  move- 
ments are  larger  than  measured  values,  indicating  the  pos- 
sibility that  the  soil  constant  determined  from  triaxial 
tests  may  not  represent  the  real  soil  strength  under  pene- 
tration (paragraphs  62  and  63). 


e. 


The  apparent  discrepancy  between  computed  and  measured  pene- 
tration resistance  suggests  that  the  usual  mathematical  pro- 
cedure, which  tends  to  magnify  the  errors  inherent  in  the 
basic  test  data,  is  not  appropriate  for  the  problem  of  soil- 
vehicle  interaction,  unless  an  objective  method  to  smooth 
the  experimental  cur”  is  developed  (paragraphs  GU  and  65). 
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Table  1 


Laboratory  Soil  Sample  Data 


Penetrometer 

Moisture 

Satura- 

Keci 

stance 

Test 

Ho. 

Type 

Speed 

cm/sec 

Content 

if 

/* 

tion 

ef 

i* 

Total 

N 

Unit 

kH/m2 

1 

5.1-cm  circular 
plate 

0.004 

49.1 

97.5 

110 

54.2 

2 

5.4-cm  circular 
cone 

0.004 

50.6 

97.5 

133 

57.6 

3 

3.6-  by  22.8-cm 
plate 

0.004 

51.7 

99.0 

350 

40.2 

It 

5.1-  by  25.4-cm 
plate 

o.oo4 

49.5 

97.0 

501 

38.6 

5 

5.1-  by  25.4-cm 
plate 

3.0 

49.7 

, 99.0 

750 

58.6 

6 

5.1-cm  circular 
plate 

0.9 

50.0 

99.0 

120 

59.  S 

7 

5.1-cm  circular 
. plate 

3.0 

49.9 

97.5 

l4o 

68.9 

8 

5.1-cm  circular 
plate 

5.6 

51.2 

97.0 

24  0 

118.5 

9 

5.4-cm  circular 
cone 

0.9 

51.0 

99.0 

l6o 

68.9 

10 

5 . 4-cm  circular 

5.6 

50.0 

97.0 

300 

129.5 

cone 
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DISPLACEMENT  BENEATH  A 5.4-CM  CIRCULAR  CONE;  PENETRATION  SPEED:  5.6  CM/SEC 


waia«MM  \m  1 : 

r ! nil  1 


lii ' IPP  *•  ■*■  ^ B C 

I . . iv" 

k ■..  : ~ 

I 

t«|»K 

gU  n 1 ' 


III  11*1.1?*!'  :i  -f  i i 

nnUlili^l  '; 

«i..«ns i»» \ i 

pHMMli  ■ ( 4'  i 

M«W9i II8t:#5wi • if i ft "i 

raMg%^sti*w!  t*  / 1 

■«&*.  v%\  ' : 

\$ . 

A ft 


*•» 

< ■ i ■ i 

■-  ifr  . /cii  1 
•Vj/h/H 


■I  -.j . >.  1.i  i> ■;; 


MMNBfJHf  I I! 

llitZ'il  j : I ' • 

■ aid  in.  i || .'  1 
■ s i m ;, 
■WKnrr.  'i  ;:  :« 

»aaais':^  -i  : i i f < 

■HMiNlI  <rt  ’l  i i 

K;aM»eiMki  i :-j  I 

SaaanMKai  zii-  n ! ; ; , 

ftll[ll^l m#i  M KM  * ■ 


•ammmmm 


AROUND  A 5. I -CM  CIRCULAR  DLATE;  PENETRATION  SPEED:  0.004  CM/SEC 
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PLATE  7 


:'JLAR  3 LATE;  PENETRATION  SPEED:  0.004  CM/SEC 
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Appendix  A:  Computer  Program 


THE  GE-4UU  SERIES  - FORTRAN  ASA  (DAPS) 
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CHOU 


f 


I 


9 

10 

11 

12 

16 

14 

15 

16 
1/ 
IB 
ly 
20 
21 
22 
26 
24 
?b 
26 
2/ 
28 
2v 

30 

31 

32 

33 

34 

35 

36 
3/ 
30 
3V 
<50 
41 
4 2 


c Y.T.CHOU  MOBILITY  RESEARCH  branch,  W.E.S. 

C PURP0SE:C0MPUTAT10N  of  strain  rate  field  in  the  soil  medium. 

DIMENSION  R(21)»Z(21)«A1(9)»A2(9)»B2(9)»A3(9),B3(9)»A4(16). 
184(16).  A(9),B(9),VRl<21).VZl(21).VR2(16.21).VZ2<16»21>.VR<16, 

1 2l),VZ(16.21).RR(16,21),CEDA(16.21).VEL(16,2l) .VRR116, 

121 >,VRZ( 16,21 ) ,VZR( 16.21 > ,VZZ<16, 21 >, 1P3<15),  ZA < 15 ) . TB < 15 > . ZB < 
lib ) » RA ( 15 ) » Y A ( 21 ) » ZC ( 15 ) • Ab ( 16 ) , Bb ( 16 ) ,IP4(15> 

COMMON  X(21),Y(21),C(16),M, I P, LP. NR , NRl , CHECK 

EQUIVALENCE  (VZ2(16,21),VZ(16.2l),VZZ(16.21)),  ( VR2 < 16 . 21 ) . VR < 16 . 
121))«  (RR(16.21).VEL(16.21).CEDA<16.21).VRZ(16.21))'.  (ZA<15). 

1TB (lb)  .ZB (15) .RA(15) ,2C(lb) ) 

READ  9000,  JJ1 
DO  3000  JJ=1,JJ1 

PRINT  1020 
1020  FORMAT  ( 1 HI , 

1 4 0Xj  48HPURP0SE  .’COMPUTATIONS  OF  SOIL  STRAIN  RATE  FIELD  //) 

READ  1 , N 1 1 , N 1 2 , N 2 A , N 2 B , N 2 c , K 1 , I P , LP , NR , NRl , RES , V , CONE  1 , CONE 2 , 

1 SHE AR1 

1 FORMAT  (10I5,f5.1.Fl0.8,3F5.1) 

READ  9000,  N2AA.N28B»NSaVe.NSaVEa»NSaVEE.N01,N02.N03»N04 .NHA.NHB, 

1 NHC.NHD  ) 

read  9000.  NH1. NH2.NH3.NH4.NHb.NH6, NUMBER. I PI . IP2 
9000  FORMAT  (1615  ) 

READ  2001,  DIST1.DIST2, AREA!, AREA2. SQUARE 
2001  FORMAT  (6F10 .5) 

IF  (CONE1-D  1001,1002.1003 

1001  IF  (CONE 2-1)  1004,1005,1005 

1004  PRINT  1010 

1010  FORMAT  (20X,40HPROJECTORi2  IN  BY  10  IN  SQUARE  PLATE  ) 

GO  TO  1014 

1005  PRINT  1011 

1011  FORMAT  (20X,45HPR0JECTORi ) .5  IN  BY  9 I N ■ SQUARE  PLATE  ) 

GO  TO  1014 

1002  PRINT  1012 

1012  FORMAT  (20X»45HPROJECTOR'2  IN  CIRCULAR  PLATE  ) 

GO  TO  1014 

1003  PRINT  1013 

1013  FORMAT  (2UX.45HPR0JECT0R!2  IN  CIRCULAR  CONE  > 

1014  CONTINUE 

PRINT  1015  , V, NUMBER 

1015  FORMAT  ( 2 0X»29HPENET RATION  SPEED ( I NCH/SEC . ) =F 1 0 . 8 , 10X , 15HTEST 


INUMBER- 1 p 
T J Nc  = SOUARE/V 
sum-o. 
do  1000 

DO  81  Ml-1'16 

DO  62  MP-1,21 

V R 2 ( M 1 , K 2 ) - 0 . 
VZ2<Ml,M2)*0. 

KR ( Ml , M2 ) = 0 , 

82  CONTINUE 
81  CONTINUE 

DO  296  M 3= 1,21 
VRl < M3 ) =0 • 


) 


s)  = 1 » 2 


AT 


•t 


THE  GE-4U0  SERIES  - FORTRAN  ASA  (DAPS) 


PAGE  # 


43 

VZ1(M3)=0, 

44 

296 

CONTINUE 

45 

IF  (J-l)  1 

46 

1030 

N1=N11 

4/ 

NSAVEBsNOl 

48 

PRINT  212 

49 

212 

FORMAT  </9X 

50 

GO  TO  1032 

51 

1031 

Nl  = Nl«? 

52 

NSAVEBsN02 

53 

PRINT  214 

54 

214 

FORMAT  </9X 

55 

1032 

CONTINUE 

56 

READ  9000, 

5/ 

READ  9000, 

58 

DO  2 

59 

IF  <ll-2> 

60 

100 

IF  (11-1) 

61 

102 

N2=N2A 

62 

N22=N2AA 

63 

GO  TO  104 

64 

103 

N2=N2B 

65 

N22=N2UB 

66 

GO  TO  104 

6/ 

101 

N^=N*C 

68 

N22=N2C 

69 

104 

Read  1044, 

70 

READ  1U44, 

71 

1044 

FORMAT 

72 

PRINT  5044, 

73 

5044 

FORMAT  </24 

74 

PRINT  1 0 4 4 , 

75 

PRINT  1044, 

76 

M = N2 

77 

IF  (11-3) 

78 

9951 

1 P= I PI 

79 

LP=I  p 

80 

GO  TO  995 

81 

9950 

IF  (Ii-4) 

82 

9953 

IP- I PI 

33 

LP=  I P 

84 

GO  T 

65 

9952 

1 P= I P2 

86 

> LP=IP 

8/ 

9954 

CONTINUE 

68 

DO  4 

89 

X ( 1 > = 7.  C 1 > 

90 

Y( I ) = R( 1 ) 

91 

4 

CONTINUE 

92 

IF  (11-5)  9 

93 

9500 

CKECK-NHC 

94 

9501 

CONt INUE 

95 

CALL  LEAST 

96 

LPP=LP4 1 

ON  THE  RIGHT  SIDE 


THE  LEFT 


< I P3  < 1 ) , !=NSAVE,NSAVEE) 
( IP4  < 1 ) . I =1.N1  ) 

11=1, N1 
100,100,101 
102.102,103 


<R< I ) , 1=1, N2) 
<zm,i=i,N2> 
(16F5.2) 

11 

HTHE  VALUES  OF  R 
(R(  I ) , I =1,N2) 
(Z(I  ), I = 1 , N 2 ) 

9950.9951,9950 


IN  I2.11HTH  COLUMN 


9952, 9953.9952 


TO  9954 


I = 1 , N2 


9500 , 9500 , 9501 


THE  GE-4U0  SERIES  - FORTRAN  ASA  (DAPS) 


97 

DO  5 

1=1. LPP 

98 

5 

Alt  l > = C( 1 ) 

99  ’ 

IF  (11-5) 

1035.1035,1036 

100 

1035 

PRINT  1037, 

(Al( I >, 1=1, LPP) 

101 

1037 

FORMAT  (10E12.5) 

102 

1036 

CONTINUE 

103 

IF  <J-1> 

7006,7006.7002 

104 

7006 

IF  (NH1) 

7002,7002,7001 

10^ 

7001 

IF  <ll-i>> 

7000.7000. 7002 

106 

7000 

PRINT  /O  05 

107 

7005 

FORMAT  (/35H 

CHECK  STATION  NUMBER 

10a 

Z2=-0.5 

1 09 

DO  / 0 o 3 

1=1,15 

110 

Z A ( I > = Z Z + 0 .5 

111 

RZ  sA1(1)+A1(2)«ZA(I)-*A1<3)*ZA<I> 

Hi? 

IF  (LP-2) 

60,60,61 

113 

61 

DO  62 

12=3. LP 

114 

IS  =12+1 

115 

tei=akis  ) 

*ZAU  >*»I2 

116 

R7=RZ+Te1 

117 

62 

continue  - 

lie 

60 

CONT  I NU[- 

119 

PRlNT  7 0 0 4 , 

ZA ( 1 ) , RZ 

120 

7004 

FORMAT  ( 2HZ 

=F5.2.40X.2HR=F10«5 

121 

ZZ  = Z A < I > 

122 

70  03 

continue 

124 

7 002 

CONTINUE 

124 

DO  6 

1=2, N2 

125 

SUU- 1-1 

126 

X( I )=SUB«T INC 

12/ 

Y ( I >=K(  I > - FI  ( 1 ) 

126 

v a < i ) - r < 1 1 

129 

x(l)=0. 

130 

Y ( 1 ) • 0 . 

131 

6 

CONTINUE 

132 

CAUL  LEAST 

133 

LPP=LP+ 1 

134 

DO  7 

1=1. LPP 

135 

7 

A 2 ( I 3 = C ( 1 ) 

136 

Y < 1 > - 0 . 

13/ 

R ( 1 ) = 0 • 

130 

DO  8 

1 = 2 , N2 

139 

SUB* I -1 

14  0 

Y ( 1 ) = S U [) « 0 , 5 

-7(1) 

141 

8 

CONTINUE 

142 

CHECK--NHC 

143 

CALL  LEAST 

144 

LPP=LP+1 

1 4 5 

DO  9 

l = 1 • LPP 

146 

9 

B 2 ( 1 ) = C ( 1 ) 

14/ 

If  < J- J > 

7015, 7015. 7010 

l4fci 

701^> 

If  < N H 2 ) 

7 0 1 0 , 7 o 1 0 . 7 0 1 1 

14  9 

7013 

If  ( 1 1 - 3 ) 

7030. 7012,7010 

150 

7012 

PRINT  7U1 9 

ONE  : 1 


AT 
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/0 1 9 


9980 


/I 


72 

/0 


7014 


7013 

7010 


292 

293 


291 


51 


i>  3 
50 
10 


11 


Vi,  OP 
Vi-03 


<U,W 


) 


) 


FORMAT  (//  4i>H  CHECK  STATION  NUMBER  TWO  S2 
PRINT  9900,  (XU),  1=1, N2) 

PRINT  990U.  < Y A < 1 ). 1 = 1#N2> 

PRINT  9980.  ( Y( I ) . I =1 • N2  ) 

FORMAT  (12E10.3) 

T T = - T INC/2. 

DO  /U 13  1=1,16 

TB( 1 >=1 T ♦T1NC/2. 

WA  = B2(1)*B2(2)oTBU>+B2(3)<»TBC1>»o2 
IF  (LP-2)  70,70,71 

DO  72  12=3, LP 

. IRR  =12*1 

TE1  = 82(IRR  ) * TB ( I ) * 0 I 2 

wa=wa+iei 

CONTINUE 

continue 

PRINT  / 0 14 , TB  ( I > , W A 
FORMAT  <2HT=F5.2,40X,2HW=F10.4 
TT  = TB ( I ) 

CONTINUE 
CONTINUE 

IF  ( N H A ) 291,291,292 

PRINT  29 3,  J,  11 

FORMAT  (//60X, 16HURAIN  WASH  AT  j=I2.5X,  6HC0LUMN  13,  8HU  THEN  «/) 
PRINT  9980 , (X(I),  1=1, N2) 

PRINT  990Q  < ( Y A ( 1 ) , ] = 1 » N2 ) 

( Y ( I ) , 1 =1 » ^2  ) 


> 


.) 


1=2, N2 


12=3, LP 


PRINT  9980, 

CONTINUE 
DO  10 

$ U B = I ” 1 

Y R 1 ( I )=A2(2)+2.«A2(3)eTJNC*SUU 

VZ1 ( I )=B2(2)+2.«B2(3)*TINC«SUB 
JF  (LP-2)  i>0  < SO  * 51 

1)0  S3 
101=12-1 

102=12*1 
T = I 2 

TE1  = T*A2(  I02)*(T1NC<*SU0)**>IQ1 

TE2  = TeB2(  IQ2)MTINC*SUB)**IU1 

VRK  I > = VR1 C I ) + T El 

VZi < I ) = V Z 1 ( I ) + TE2 

CONTINUE 

CONTINUE 

CONTINUE 

VR  1 ( 1 ) = 0 , 

VZ1(1)=0- 

DO  11  I =1 , N2 

x < i ) = z ( i ) 

Y ( I ) = V H 1 ( I ) 

Y A ( 1 ) = Y ( I ) 

CONTINUE 

IF  (11-1)  9502, 9502,9503 

CHECK=NHC 

CONTINUE 


i 
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702S. 7U2S. 7020 
7020. 7(120. 7021 
7020. 7U22. 7020 


M = N2 

CALL  LEAST 
LPP=LP+1 

DO  12  . 1=1, LPP 

A3( I )=C( I ) 

IF  (11-1)  6000.6000,6001 

PRINT  10.57,  (A3(  I ),  1=1, UPP) 

CONTINUE 

DO  13  1=1, N2 

Y< 1)=V^1< 1) 

CALL  LEAST 
LPP=LPf 1 

DO  14  1=1, LPP 

B3( I )=C<  I ) 

IF  (J-l>  ‘ 702S.702b.7020 

IF  (NH3)  7020.7020,7021 

IF  ( H-3)  . 7020. 7U22. 7020 

PRINT  /029 

fOHMAT  (//  4 OH  CHECK  STATION  NUMBER  /THREE : 3 (VR  V Z ) 

PRINT  9980,  ( X ( 1 ) . 1=1,  N2) 

PRINT  WOO,-  ( YA(  1 ) . 1 = 1, N2) 

PRINT  9980.  (Y(I>, 1=1, N2> 

ZZ  = -1. 

DO  7023  1 = 1,10 

7B ( I >=ZZ+1. 

VV=B3(l)+B3(2)ftZB( I ) +B3 ( 3 > *Z0 ( 1 )««2 
IF  ( LP-2 ) 86,86,87 

DO  88  12=3. LP 

101=12+1 

TE1  = U3(  I0JL)*ZB(  I )*«I2 

VV=VV+ 1 E 1 

CONTINUE 

CONTINUE 

PRINT  7024,  ZB  < I > . V V 

FORMAT  (2HZ  = FS.2,4GX,24HVEL0(Jl  TY  IN  Z D I RECT  1 ON  = F 14 . 8 ) 

ZZ=ZH( 1 ) 

CONTINUE 

CONTINUE 

IF  (NHB)  294,294,708b 
PRINT  7086  . j,n 

FORMAT  (//60X.16HBRAIN  HASH  AT  J=I2.bX,  6HC0LUMN  I3.9HVR  AND  VZ/) 
PRINT  9980,  (X(I).  1=1, N2) 

PRINT  9980  , ( YA ( I ) » 1 =1 » N2 ) 

PRINT  9980.  ( Y < I ), 1=1, N?) 

CONTINUE 

DO  li>  12  = N$AVE . N22 

SUB =12-1 

VR2  ( 1 1 , 1 2 ) = A3  ( 1 ) +A3  ( 2 ) * SU8»0  . S + A3  ( 3 ) « ( ( SU8‘>0  . b ) <>«2 ) 

VZ2  ( 1 1 » I2)=B3(l>+B3<2>*SU3«0.S  + iT3t3>«(  ( SUB+' 0 • S ) « *?  ) 

RR(  11 . 12  )=A1  ( 1 > + A1  <2)  tSl'B  = 0 ,b  + Al  (3)*  ( ( SUB«  0 . ) ««2  ) 

IF  ' Lp-2  ) 90,90,91 

DO  92  I = 3 , L P 

101=1+1 

T E 1 = A3 ( I01)c ( SU8°C • b ) « « I 


1 2=NSA VE . N22 


Sfc,  •.< 
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TE2=B3 ( I 01 ) • ( SUB«  0 . b ) 0*  I 

260 

TE3  = A1(  I Q1 ) • ( Sl)B<*0 , 5 ) **  I 

261 

VR2  < 11. I 2 » = VR2  t 11. 1 2 1 ♦ TE1 

262 

VZ2<  U,  12 ) =VZ2 ( 11, 1 2 ) *TE2 

264 

RR(I1, l2)=RR(ll. I 2 ) ♦ TE3 

264 

92 

CONTINUE 

26b 

90 

continue 

266 

IF  <KR(1. 121-0.0000001)  301,301,302 

26/ 

301 

RR ( 1 . 1 2 ) =0 , 

26B 

302 

CONTINUE 

26V 

15 

continue 

2/0 

UO  4000  1=1. N22 

2/1 

VR2 ( 1 » 1 ) = VR2 ( 2 « I ) 

272 

VZ2 ( 1 , I ) = VZ2 ( 2 > I ) 

2/4 

4000 

CONTINUE 

2/4 

IF  (J-l)  7036.7046,7040 

27b 

/036 

IF  ( NH4  ) 7030.7040.7041 

276 

/041 

IF  (11-4)  7032, 7042, 7040 

2// 

/042 

PRINT  7039,  11 

2/6 

7049 

FORMAT  (//  52H  CHECK  STATION  NUMBER  FOUR: 

1 

1 5 , 2 0 H T M COLUMN 

279 

PRINT  7U35,  ( RR ( 1 1 . I 2 ) , I 2-NSAVE . N22  ) 

280 

PRINT  /04b,  <VR2( 11, 12), I2=NSAVE,N22) 

281 

PRINT  /03b,  <VZ2( 11, ]2>, 1 2=NSA VE , N22 ) 

282 

/045 

FORMAT  < 1 2 E 1 0 • 3 ) , 

284 

7030 

continue 

284 

2 

CONTINUE 

28b 

PRINT  /0/3 

286 

70/4 

FORMAT  ( 4 0 X , 3 0 H LOOP  2 PRINT  OUT  COMPLETED 

28/ 

N2=NS A VEE 

288 

no  20  12=NSAVE,N2 

28V 

IP=IP4( 12) 

290 

• 

LP  = IP 

291 

IF  (I2-N2AA)  200,200,201 

292 

200 

MA  = 1 

294 

GO  TO  202 

29  4 ■ 

201 

IF  (I2-N2BB)  203, 20  31,  204 

29b 

203 

N A =2 

296 

GO  to  202 

29/ 

204 

NA=4 

298 

202 

CONTINUE 

299 

M=N1-NA+1 

400 

DO  21  1 1 = N A » N 1 

401 

I = U-(NA-1  ) 

402 

X(  l ) - R R ( 11.12) 

404 

Yd  )=VR2(  11,12) 

404 

Y A ( I ) = Y ( 1 ) 

40b 

21 

CONT 1 NUE 

406 

IF  (I2-N2AA)  9 b 0 b , 9 b 0 b , 9 b 0 4 

40/ 

9b  Ob 

check=nhc 

408 

Vb  0 4 

cONti Nut 

409 

CALL  LI  A ST 

410 

l.  PP-LP*1 

411 

DO  22  I =1 , LPP 

J 


) 
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I1=NA,N1 


A4 ( I ) =C ( l ) . 

If  ( I 2-NiJ  A A ) 6010.6011.6010 

PRINT  1037,  <A4<1), 1=1. LPP) 

CONTINUE 

DO  23  I1=NA,N1 

I = U-(NA-1) 

Y ( I )=VZ2(  11.12) 

CHECK=NHC 
CALL  LEAST 

lpp=lp+i 

DO  24  1=1* LPP 

B4( I )=C( 1 ) 

IF  (J-l)  7043,7043.7040 

IF  <NH5>  7040.?040.7041 

IF  (I2-N2AA)  305.305.  7040 


I =1 • LPP 


/047 

9980, 


( RR ( 1 1 » I 2 ) » 11  = 1, 

« X ( I ) . 1=1, Nl) 

(YA(  1 ), 1 =1,N1  ) 

( Y ( I ) , I = 1 . N 1 ) 


1=1.16 


•11  = 3- LP 


PRINT  7045, IP 

FORMAT  ( / /21HCHECK  S1ATION  NO. 5=  13) 

PRINT  7046,  12 

FORMAT  ( 5 0 H VELOCITY  IN  l.  DIRECTION  ALONG  HORIZONTAL  LINE-  12=12 
PRINT  7047-,  ( VZ2 ( 1 1 , I 2 ) , 11  = 1. Nl  ) 

FORMAT  (12E10.3) 

PRINT  /048  > 

FORMAT  ( 5 0 H THE  CORRESPONDING  RR , THEN  CORREPONDING  VR  AND  VZ 
PRINT  /\)4!  , ( R R ( 1 1 » I 2 ) , 11  = 1, Nl  ) 

PRINT  9980,  ( X ( I ) , 1 =1 , N)  ) 

PRINT  WOO,  ( YA  ( 1 ) , 1 =1,  Nl  > 

print  9900.  < y< i > . i =1 ,ni ) 

RAA=-0 .25 

00  7050  1=1,16 

RA  ( I ) = F(AA*0 , 25 

VRA=A4(l)+A4(2)6RA(!)+A4T3)«RA(I)e*2 
VZA=B4(l)*B4(2)*RA<I)*U4(3)eRA(I)*«2 
IF  (LP-2)  95.95,96 

DO  'I  1 = 3 > L P 

101=11*1 

T El = A 4 ( I Q 1 ) e R A ( 1 ) *» I 1 
T E 2 = H 4 ( I Q 1 ) o R A ( I ) »«  1 1 
VRA  = VRA*  T E 1 
VZA  = VZA  . IE2 
CONT INUE 
CONT I NUE 

PRINT  /051,  RAC  I J.VRA.VZA 

FORMAT  (2HR  = F5.2»20X»4HVR=  E12 . 5 . 20X  , 4HVZ  = 

R A A = R A ( I ) 

CONT I NOE 
CONT INUE 

DO  25  1 1 = N A » N 1 

SUU= 1 1-1 

VRC 1 1 . I2)  = A4 ( l ) * A 4 ( 2 ) « SUB » 0 . 5 * A 4 ( 3 ) o ( SU 9 » 0 
V Z ( 1 1 , I 2 ) = B 4 < 1 ) +B4(2)»SUB*0 .5*84 ( 3)<(SUC* 0 
IF  (l.P-2)  120.120.121, 

DO  J22  1 = 3 , L P 

101=  1*1 

T f:  1 = A 4 ( 101)«(SUB«0.5)**I 


E12.5) 


5 ) « » 2 
5 > * * 2 


A 7 


uij  VfpngpmBi 
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i 

I 

i 

I 


>566 
36  / 

3 68 
469 
4/0 
4/1 
3/2 
4/4 

4 7 A 

4 7 1> 
4/6 
4// 
4 78 
479 
48  0 
481 
382 
464 
48  4 
38V 
486 
48/ 

488 

489 

490 

491 

492 
4 94 
494 
49V 
496 
49/ 

498 

499 

400 

401 
4 0 <? 
404 
404 
40V 
406 
40/ 

408 

409 

410 

411 

412 


414 

414 

41V 


122 

120 

2 V 

26 

9002 

9996 

9997 

27 

28 
29 


TE2=B4 ( JQ1)»(SUB»0.V)*«I 
VR( !1» I 2 ) = VR ( 11» I 2 ) ♦ T L 1 
VZ<  11.  1 2 ) = VZ ( 11.  1 2 ) + I F:  2 ' 

CONTINUE 

CONTINUE 

VELdl.  12)=S0RI  (VR(11.  I2)»»2+VZ<  11.  I2)**2) 

CONTINUE 

PRINT  26,  12 

FORMAT  (/  48H  THE  COMPUTED  VELOCITIES  AT  GRIU  POINTS  ALONG 

1 I 3 • 44h  Th  ROW ( COUNTED  FROM  THE  BOTTOM  OF  THE  MOLD)  ) 


111 


PRINT  9002,  (VELdl,  12),  11  = 1, 
FORMAT  (12E10.3) 

PRINT  9998 
FORMAT  (lOX.VHVZ! 

PRINT  9002,  <VZd 

PRINT  9997 
FORMAT  (1CX.,VHVR: 

PRINT  9 U 0 2 , (VRd 

no  2/ 


N1 ) 


I 2) , I=l,Nl 


I 2 ) , I = 1 , N 1 
ii=na«ni 


) 


CEDA(  II,  I2)=ATAN(VZdl,  I2)-/VR(  II.  12)  )«57.2 
PRINT  28 

FORMAT  (44HTHE  DEGREE  OF  ANGLE  AT  SAME  ROW 
PRINT  29,  ( C E D A ( II,  I2)».I1  = 1>N1) 

FORMAT  (16F7.1) 

DO  4 0 1 1 £ N A » N 1 

SUB=I1-1 

VRRt 11’ 12)=A4(2)+2i«A4(3)»SUBo0"V 
VZR( II* I2)=B4<2)+2.*B4(3)«SUD*0.5 


112 

110 

40 

221 

222 

223 

224 
220 


IF  ( L P - 2 > 110,110,111 

DO  112  1=3, 

101=1+1 
102=1-1 
T = I 

TE1  = T*A4<  J01)<MSUB*0-V)®°J02 
TE2  = UB4<  I 01  )»  (SUB  *0.5)  **102 
VRK ( 11. I 2 ) = V R R < II, 12 ) +TE1 
VZRt  II . I 2 ) = VZR ( II. 12J  + TE2 

continue 

CONTINUE 

CONTINUE 

PRINT  222,  1 2 

FORMAT  ( 4 OH  THE  VALUES  OF 
PRINT  223,  ( VRR (11,12), 

FORMAT  (12E10.3) 

PRINT  224 

FORMAT  ( 20H 

PRINT  223,  ( VZR ( 1 1 . I 2 ) 

CONTINUE 


LP 


VRR  A 1 GRID 
1 1 = 1 < N 1 > 


VZR 

11  = 1.  ND 


POINTS  ALONG  I3.6HTH  ROW 


20  CON T I NUE 


PRINT  7071. 

7 0/1  FORMAT  ( / / 7 1 0 X 1 


10HLO0P  20  CK 


///) 


AP 


W 


THfc  GE-4UU  SERIES  - FORTRAN  As*  (OARS) 

DO  32  11=1, N1 

IP=1P4(U) 

LP=  IP 

IF  < 1 1 -i? ) 230.23u.231 

230  IF  1)1-1)  232,232.233 

232  N2=N2AA-N03 
GO  TO  234 

233  N2=N2BU-N04 
GO  TO  234 

231  N2=N^C 

234  m=N2-NSAVE+1 

DO  33  l2  = fJSAVb,N2 

SUR=I2-1 

16=l2-(NSAVc-l ) 

X<  I 6 ) =SUB»0  • 5 
Y(  I6)=VR( II, 12) 

33  CONTINUE 

IF  ( 1 1-3)  9507,95)07,9506 

1*507  C'Hfc CK  =NHC 
9t>06  CONI  INUE 

CAUL  LEAST 
LPP=LP+1 


DO  34 

34  A5<I  )sC(l ) 

IF  ( 1 1 - 3 ) 602(1.. 

&021  PRINT  1 0 3 ? , 

6020  CONTINUE 

DO  35) 

1 6= I 2- ( NS AVE - 1 ) 

Y(  I6)=VZ( 11 < 12) 

35  CONTINUE 

CAUL  LEAST 

lpp=lp+i 

DO  36 

36  B5( I >=C( 1 ) 

IF  ( 1 1-NSAVEO ) 
3011  IF  (J-l)  7q6 

/063  IF  ( N U 6 ) 70 


1=1, LPP 


6020.6021.6021 


( A5(  i > 


1 = 1 .LPP) 


7 0 6 3 , 7 o 6 3 , 7060 
7060,7060,7061 


I2=NSAVE.N2 


1=1 < LPP 

3011.301). 7060 


7061  PRINT 
7064  FORMAT 


7064  , 1 1 , IP 
(5011 


ALONG  COLUMN 


13,12 

PRINT  /062  , < V R ( 1 1 , I 2 ) , 12  = 1. N2C  ) 

PRINT  7U62  , ( VZ ( 1 1 , J 2 ) , 12  = 1. N2C.  ) 

7 () 6 2 FORMAT  (10E10.3  ) 

Z Z = - 0 , 7 1> 

DO  30U6  1=1.12 

7C( I >=ZZ*0 . 75 

RZl=A5(l)+A5(2)*2C(l)+A5(3)*ZC<l)e«2 

RZ2r£35(l  )*D5>(2)*ZC<  I ) +D5  ( 3 )<^C  ( ! )®»2 


12=1. N2C 
I 2 = 1 « N 2 C 


IF 

3000  DO 


IF  (LP-2  > 
DO  3009 
IS=I2*1 
T E 1 = A 5>  ( IS) 

T E 2 = B 5 ( IS  ) 


3007,3007,3006 
12  = 3, LP 


«7C<  I >«ol2 
« Z C ( I ) « * I 2 


A9 
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4 69 

R21=RZ1*TE1 

4/0 

RZ^sRZ«?  + TE2 

471 

3009 

CONTINUE 

4 72 

3007 

CONTINUE 

476 

PRINT  3010,  ZC( 1 ) ,RZ1,RZ2 

474 

3010 

FORMAT  (E12.2,10X.E12.2,10X>E12.5) 

47t> 

ZZ  = ZC(  1 ) 

476 

3006 

CONTINUE 

47  7 

7060 

CONTINUE 

470 

DO  67  I2=NSAVE,N2 

4/9 

SUti-lP'l 

480 

VRZ<  11 . I2)=A5<2>+2.  *>A2(3)»SUB»0 .5 

481 

VZX< ll» I2>=B2(2)*2.«82(3>*SU3e0.2 

482 

II  (LP-2)  1i5,115,h6 

4 83 

116 

DO  117  1=3,  l.P 

484 

101=1-1 

4 8i> 

j 0^= I + 1 

406 

T = 1 

40/ 

Tfcl=T«A2(lO?)*( SUB* 0 .2 )«« 1 Ql 

488 

TE2  = T»D2(  102)*  ( SUR*  0 . i> ) *«  1 Ql 

489 

VRZ  ( 11.  I 2 ) = V f?  ? ( 11,  I 2 ) + T 1 1 

490 

VZZ( 11 . I 2 ) = V Z Z ( 11 , I 2 > + T E 2 

491 

117 

CONTINUE 

492 

112 

continue 

49  6 

37 

continue 

494 

PRINT  6030,  11 

492 

6030 

FORMAT  (26H  VRZ  AT  GRID  POINTS  ALONG 

4 96 

PRINT  223,  ( VRZ ( 11 > 12) > 1 2 = NS  AVfc • 

49/ 

PRINT  6031 

498 

6031 

FORMAT  ( 20H  VZZ 

499 

PRINT  2/3 , ( V Z Z ( 11,  12),  1 2 = N S A V E , 

200 

32 

CONt I NUE 

201 

PR  1 NT  70 72 

20)? 

70/2 

FORMAT  ( 60X, 1 DHL OOP  32  OK 

203 

S U M 2 = 0 • 

204 

DO  40  U = 1,NSAVEB 

202 

IF  Ul-2)  240.240,241 

2y6 

24o 

IF  ( 11-  1 ) 242,242.243 

20/ 

242 

N 2 = N 2 A A - N 0 3 

208 

ARE  A = ARC  A 1 

209 

D I S T = D I S T 1 

210 

GO  TO  244 

211 

243 

N 2 = N 2 B 8 - NO  4 

21? 

SUB= I 1-1 

21  3 

D1ST  = 0 ,t>«SUB  + D15T2 

214 

AF!E  A = ARE  A2 

212 

CC  TO  244 

216 

24) 

N2=NSAVEt 

21/ 

SU«s I 1-1 

2 1 8 

DIST  = 0 •S»*,SUS*D1ST2 

219 

A R E A : A R [:  A 2 

22  0 

24  4 

(.O',  T I NUt 

2?  V 

PRINT  130,  1 1 , N S A V F A 

2?? 

130 

f ORMAT  ( / 4 D X , 24riTi,i:  STRtSS  COMMON;  Mf 

13.11HTH  COLUMNS 


I 3 » 9HT H COLUMN  7H:N$AY>  = 


J 
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1 13  , 

523  PRINT  131 

524  131  FORMAT  <5X,  95HSR  SC  SZ 

1 SRZ  J2  /> 

t»2t»  SUMl  = 0. 

526  IF  ( K 1 — 1 ) 261,261.260 

52 / 260  DO  41  I2  = NSAVEA,N2 

52b  IF  < 1 2-N2 ) 3001,3002,3002 

529  3002  AREA=AREA/2. 

530  3001  CONTINUE 

531  FR  = VRR<  U,  12) 

532  fcC=VR( II' 12) /D 1ST 

533  EZ^VZZt II. 12) 

534  GAMA=VKZ( 11. 12)+VZR( 11. 12) 

53&  VA1=ER**24EC*»24EZ*»2 

536  vA2=l . ♦ (GAMA»«2/( vA1iGAMA»«2) ) 

53/  VA3rSQRT ( ( 0 ,5*VAl+0 .25*GAMA*«2)/< VA1**2 ) ) 

‘->36  VA;>=  (0  .5*VAl+0  .25«GAMA«»?)/(  VAl*«2) 

539  \'A4  = 2 .*3. 14159*DI  ST«  AREA*  ( VA2/VA5)*VA3 

540  SUM1=SUM1+VA4 

541  SR  = 1./SQRT(0.5«(1 . 4<EC*e2/EK«»2)4<EZ»*2/ERa*2) ) 4 0 . 25« C G AM A« * 2/ERa * 

12  5)  > 

542  SC  = 1 . / SORT  <0 . 5*  ( <E:R»*2/ECa*2  )41 . 4 (EZe*2/EC««2)  )+0 . 25*  (GAMAa^/ECee 
12  > > 

543  SZ  = 1.  ZSQRTlO  ,5«(  (eR««2/EZ**2  )4  <ECa*2/EZ»«2>4l . )■•  0 .25*  (GAMAa*2/EZ«» 
1?  ) ) 

544  SRZ  = 1 . / S Q R T ( 2 . 0 * V A 1 / ( G A H A * « 2 ) 4 1 .00) 

5 45  WU2=  0 . 5* (SR a 5 24  SC* «24S2**2) 4 SRZ* *2 

546  IF  ( NHL) ) 41,41,200? 

54/  2002  PRINT  132.  SR . SC . SZ , SRZ , HJ2 

540  132  FORMAT  < E12.5.5X,  E 12 . 5 . 1 OX , E 12 . 5 . 7X , E 1 2 . 5 , lflX , E 1 2 . 5 ) 

549  41  CONTINUE 

550  GO  TO  2/5 

551  261  00  42  1 2 = NS AVfc  A , N2 

552  IF  ( 1 2 ~ N 2 ) 3.  3*3004  *3004 

553  3004  ARE A= AREA/2 . 

554  3003  CONTINUE 

555  ER=VRR« 1 1,12) 

556  EZ=VZ?( 11* 12) 

55/  GAMA=VRZ ( ! 1 . I?)4 VZR< 1 1 . 12 ) 

55b  VA1=ER«»24EZ**2 

559  V A *^  = 1 . 4 ( G AM  A * » ’2/  ( V A 1 + G AM  A a * 2 ) ) 

560  V A 3 = S 0 R T ( ( 0 , 5 * V A 1 4 0 . 2 5 a G A M A a * 2 ) / ( V A 1 a a 2 ) ) 

561  VA5  - (0 .5*VA1  +0 .25a  GAMA*  *2  ) / ( V A 1 a * 2 ) 

562  V A 4 - A R E A * { V A 2 / V A 5 ) a V A 3 

563  SUM1  = SUMl 4 V A 4 

564  SR  * 1 . / S 0 R T ( 0 , 5 * ( 1 . 4 ( E Z * * 2 / E R * * 2 ) ) 4 0 . 2 3 * ( G A N A * * ? / E R * * 2 ) > 

565  SC  *0. 

566  SZ  a 1 . / S 0 R T ( P , 5 * ( ( ! R a z / E Z * * 2 ) 4 1 . )40,?5»(GAMA4*?/LZ**2>) 

56  7 SR  Z a 1 . / SORT ( 2 . 0 » V A1 / ( U AM  A ^ *2  ) 4 1 . 0 0 ) 

566.  KJ2  sO . 5» ( SR  «*24S7  **2)4SR2  a«? 

569  IF  (NHL’)  4 7,42.2003 

5 Zy  2003  PRINT  2/3,  SR, SC » SZ » SRZ . W J2 

5/1  2/3  FORMAT  ( E 1 ? • 5 , 5 X , El  ? . 5 , 1 OX  , E 12 . 5 . 7 X > E ! ? . 5 , 1 6 X , E- 1 ? . 5 ) 
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572  42  CONTINUE 

57 3 2/5  PRINT  43,  11, SUM1 

>74  4,5  FORMAT  ( / 30X,  10HTHE  SUM  OF  I3,7HC0LUMN=E12.5 

575  SUM2=SUM2+SUM1 

576  40  CONTINUE 

577  IF  (NMD)  2004.2004,2005 

578  2005  PRINT  45,  SUM2 

579  45  FORMAT  (32HTHE  TOTAL  SUM  OF  THIS  HALF  M0LD=E12.5 

580  2004  SUM=SUK+SUM2 

581  1000  CONTINUE 

582  ENL=REScV 

583  PRINT  46*  RES'SUM'EnL 

584  46  FORMAT  (20X.30HTHE  MEASURED  RESISTANCE  F0RCE=F10.2  / 

1 20X , 30HTHE  COMPUTED  ENERGY  =El2.5  / 

1 20X  *26H  COMPUTED  R«V  = E12.5/ 

585  3000  CONTINUE  • 

586  CALL  EXIT 

58/  END 


at: 


///) 
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1 

suhrouiine  least 

2 

COMMON  X(21)»Y(21),C(16).M 

3 

P I MENS  1 ON  P IM1 ( 50 ) • P 1 < bO  ) i 
1 PtRH ( 50 ) , T ( 3 ) 

A 

1 (1 )=2.BbubGS2241 

5 

1 (2)  = -2.0/514  80497E-03 

6 

1 <3)=-2.2b20J  49609E-01 

/ 

7 b R = 0 • 

8 

ONt  = l . 

V 

IF  (LP-IF*)  104,106.106 

10 

104 

l P=  I P 

11 

1 06 

N = L P + 1 

12 

K = 1 P + 1 

13 

120 

PO  130  J=1,M 

14 

P 1 Ml ( I l^ZbR 

lb 

PI ( 1 )=ONb 

16 

130 

CONTINUE 

1 / 

[’0  14  0 1=1,16 

IB 

S<  1 )=ZER 

19 

B L » A ( I >--ZbR 

20 

140 

ALpHA( 1 ) - Z E R 

21 

A1=ZLR 

22 

pl=ZLR 

23 

R1I=M 

24 

1=1 

2b 

16  0 

Rl=ZtR 

26 

D0‘  1/0  L = 1,M 

2/ 

1 /o 

RIsWJ+Y(t.)*Pl(L) 

26 

S(1  ) = W]/W1 1 

29 

IF  Cl-N)  180,302,302 

30 

1 GO 

A 1 = Z - R 

31 

DO  190  L=1,M 

32 

190 

A1=A1+X(L  >«PI  <L  )**2 

33 

Al= Al/W 1 1 

34 

ALPHA! 1 + 1 ) = A 1 

3b 

R I = Z b R 

36 

DO  210  L = 1 , M 

3/ 

P I PI  = ( X ( L ) " A 1 ) *P  i ( L ) "HI  * 1’  I M: 

3B 

PJMKI.  ) = B 1 (L) 

39 

PI  (L)=P1P1 

40 

N 1 =W  1 *P I PI »P 1 PI 

41 

211; 

CONI  1 NUB 

42 

D 3 W 1 / W J 1 

4 3 

BE  1 A < I + 2 ) = R 1 

44 

W I 1 s R 1 

4b 

1 = 10 

4 6 

00  TO  160 

4/ 

302 

DO  310  L = 2 » 1 6 

4 B 

310 

G < L ) = Z t R 

49 

G < 1 ) --  Of  . B 

DO  3bC  J = 1,K 

S 1 = Z r H 

|:0  320  L = l,p 

G(L)=G<L  ) - ALPHA  (L)>G<1  - 1 > -UL 7 A (L > °G< l - ?> 
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54 

S1=S1*S(L)*G(U 

55 

320 

CONTINUE 

56 

C( J)=S1 

57 

L=K 

58 

DO  340  12  = 2*  K 

59 

G(L)=C(L-1) 

60 

340 

L=L-1 

61 

350 

G(1  )=ZER 

62 

If  ( NR1 ) 404.405.404 

6 A 

404 

PRINT  V05,  ( J-l ,C( J) , J=1,K) 

64 

905 

fORKAI  (8E12.5) 

65 

405 

CONTINUE 

66 

IF  (CHECK)  1050.1050,1051 

6 7 

1053 

PRINT  1035 

68 

1035 

f OWNA  T ( 26HC0EF  F I C I ENT  S ARE  *** 

6^ 

LPP  = 1.P  + 1 

7 0 

PRlNl  1036,  (C(I).1=1,LPP) 

71 

1036 

FORMAT  (10E12.5) 

) 


| 

I 

! 

i 

I 


i 


| 

I 
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